** Make figure 3, showing effect of an additional $1000 of EITC exposure
** on the probability of being in each of five categories of self-reported
** health status
****************************************************************************


gen Health1 = (ehltstat==1)
replace Health1 = . if ehltstat==0
gen Health2 = (ehltstat==2)
replace Health2 = . if ehltstat==0
gen Health3 = (ehltstat==3)
replace Health3 = . if ehltstat==0
gen Health4 = (ehltstat==4)
replace Health4 = . if ehltstat==0
gen Health5 = (ehltstat==5)
replace Health5 = . if ehltstat==0

areg Health1 EITC_v3 tage AgeSq Female Married White HSGrad SomeCollege CollegeGrad  unemp_rate r_MaxWelf state_rate_wages minwage gdp any_w_a_extended i.rhcalyr i.N_PotentialDep if (NoParsWithCollege==1 & MissingHealthInfo==0), absorb(tfipsst) cluster(tfipsst)
estimates store Health1

areg Health2 EITC_v3 tage AgeSq Female Married White HSGrad SomeCollege CollegeGrad  unemp_rate r_MaxWelf state_rate_wages minwage gdp any_w_a_extended i.rhcalyr i.N_PotentialDep if (NoParsWithCollege==1 & MissingHealthInfo==0), absorb(tfipsst) cluster(tfipsst)
estimates store Health2

areg Health3 EITC_v3 tage AgeSq Female Married White HSGrad SomeCollege CollegeGrad  unemp_rate r_MaxWelf state_rate_wages minwage gdp any_w_a_extended i.rhcalyr i.N_PotentialDep if (NoParsWithCollege==1 & MissingHealthInfo==0), absorb(tfipsst) cluster(tfipsst)
estimates store Health3

areg Health4 EITC_v3 tage AgeSq Female Married White HSGrad SomeCollege CollegeGrad  unemp_rate r_MaxWelf state_rate_wages minwage gdp any_w_a_extended i.rhcalyr i.N_PotentialDep if (NoParsWithCollege==1 & MissingHealthInfo==0), absorb(tfipsst) cluster(tfipsst)
estimates store Health4

areg Health5 EITC_v3 tage AgeSq Female Married White HSGrad SomeCollege CollegeGrad  unemp_rate r_MaxWelf state_rate_wages minwage gdp any_w_a_extended i.rhcalyr i.N_PotentialDep if (NoParsWithCollege==1 & MissingHealthInfo==0), absorb(tfipsst) cluster(tfipsst)
estimates store Health5

coefplot (Health1)  (Health2)  (Health3)  (Health4)  (Health5), keep(EITC_v3) vertical legend(off)  ytitle("EITC Effect") msymbol(o) lcolor(black) coeflabels(EITC_v3="Self-Reported Health Category") scheme(s1mono) yline(0)
graph export Health_Category_Figure.eps